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1 Introduction 



Recently [1] we have suggested some kind of a coalescence model [2] for the description 
of the hadronization from the quark-gluon plasma (QGP) phase of QCD considered as a 
thermalized quark-gluon system [3] at high densities and temperature in which quarks, 
antiquarks and gluons being at the deconfined phase collide frequently each other. There 
is a belief [4] that the QGP phase of the quark-gluon system can be realized in ultra- 
relativistic heavy-ion collision (£' cms /nucleon 1 GeV) experiments. 

At very high energies of heavy-ion collisions the quark-gluon system is composed 
from highly relativistic and very dense quarks, antiquarks and gluons. By virtue of the 
asymptotic freedom the particles are almost at liberty and due to high density collide 
themselves frequently that leads to an equilibrium state. If to consider such a state as a 
thermalized QGP phase of QCD, the probabilities of light massless quarks n q (p ) and light 
massless antiquarks riq(p), where q = u or d, to have a momentum p at a temperature T, 
can be described by the Fermi-Dirac distribution functions [3,5]: 

= e -u(T) lp/T ~ I ' n 'W = e v(T)+ P /T + 1 > (L1) 

where a temperature T is measured in MeV, v(T) = /x(T)/T, /x(T) is a chemical potential 
of the light massless quarks q = u,d, depending on a temperature T [5]. A chemical 
potential of light antiquarks amounts to — /i(T). A positively defined fJ.(T) provides an 
abundance of light quarks with respect to light antiquarks for a thermalized state [1,4]. 
A chemical potential fi(T) is a phenomenological parameter of the approach which we 
would fix below [1]. 

The probability for gluons to have a momentum p at a temperature T is given by the 
Bose-Einstein distribution function 

n ' (p) = jihi- (L2) 

Since a strangeness of the colliding heavy-ions amounts to zero, the densities of strange 
quarks and antiquarks should be equal. The former implies a zero-value of a chemical 
potential /i s — /ig — 0. In this case the probabilities of strange quarks and antiquarks can 
be given by 

n 8 (p) = n s (p) = 1 , (1.3) 

e y/p 2 + ml/T + x 

where m s = 135 MeV [6] is the mass of the strange quark and antiquark. The value of 
the current s-quark mass m s = 135 MeV has been successfully applied to the calculation 
of chiral corrections to the amplitudes of low-energy interactions, form factors and mass 
spectra of low-lying hadrons [7] and charmed heavy-light mesons [8]. Unlike the massless 
antiquarks u and d for which the suppression is caused by a chemical potential /i(T), the 
strange quarks and antiquarks are suppressed by virtue of the non-zero mass m s . 

In Ref. [1] we have supposed that a chemical potential /x(T), a phenomenological pa- 
rameter of the description of the QGP state as a thermalized quark-gluon system at a tem- 
perature T, is an intrinsic characteristic of a thermalized quark-gluon system. Thereby, 
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if the QGP is an excited state of the QCD vacuum, so a chemical potential should exist 
not only for ultra-relativistic heavy-ion collisions. Quark distribution functions of a ther- 
malized quark-gluon system at a temperature T should be characterized by a chemical 
potential fi(T) for any external state and any external conditions. Since any state of 
a thermalized system is closely related to external conditions, in order to obtain fi(T) 
we need only to specify the external conditions of a thermalized quark-gluon system the 
convenient for the determination of /x(T). 

Indeed, it is well-known [5] that the Helmholtz free energy F(T, V, N) defining the 
partition function Z(T, V, N), F(T, V, N) = —T £nZ(T, V, N) which plays a central role 
in studying thermalized systems, is nothing more than a work for an isothermic process. 
Therefore, by producing external conditions keeping T = const and measuring a work 
one can get a full information about the Helmholtz free energy F exp (T, V, N). Then, in 
terms of this Helmholtz free energy F exp (T, V,N) one can obtain the partition function 
Z Pm (T, V, N) which can be applied to the description of the thermalized system at any 

Following this idea we have fixed the chemical potential fi(T) in the form [1]: 
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where /io = /x(0) = 250 MeV is a chemical potential at zero temperature [1]. The In- 
dependence of a chemical potential given by Eq.( |1.5|) differs by a factor 1/4 from the 
low-temperature behaviour of a chemical potential of a thermalized electron gas [9] . The 
former is caused by the contribution of antiquarks. 

In the high-temperature limit T — > oo a chemical potential /x(T) defined by Eq.(L4) 
drops like T~ 2 : 



1.6) 



A chemical potential drops very swiftly when a temperature increases. Indeed, at T = 
160 MeV we obtain fi(T) ~ £*o/4, while at T = /x a value of a chemical potential makes 
up about tenth part of /io, i.e. /x(T) ~ /io/10. This implies that at very high temperatures 
the function v(T) = fi(T)/T becomes small and the contribution of a chemical potential 
of light quarks and antiquarks can be taken into account perturbatively. This assumes in 
particular that at temperatures T > /x = 250 MeV the number of light antiquarks will 
not be suppressed by a chemical potential relative to the number of light quarks. 

In our approach the multiplicities of hadron production we define in terms of quark and 
anti-quark distributions functions in a way similar to a simple coalescence model [2] but 
for correlated quarks and anti-quarks. Indeed, in a coalescence model quarks and anti- 
quarks are uncorrelated [2]. This allows to introduce separately the number of light quarks 
q and light anti-quarks q and the number of strange quarks s and strange anti-quarks s 
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[2]. The subsequent calculation of multiplicities of hadrons in a simple coalescence model 
resembles a quark counting. In fact, multiplicities of hadrons are proportional to products 
of the number of quarks (q, s) and anti-quarks (q, s) in accord the naive quark structure 
of hadrons. Since quarks and anti-quarks do not correlate, so that the multiplicities of 
hadrons turn out to be independent on the momenta of hadrons []. 

In our approach the multiplicities of hadrons produced from the QGP phase are de- 
scribed by momentum integrals of quark and anti-quark distribution functions. Unlike a 
coalescence model these integrals depend explicitly on the momenta of hadrons, temper- 
ature T and chemical potential /x(T). 

For example, the multiplicities of the production of the mesons, NK±(q,T) } and 
^ mesons, N w ±(q,T), are defined in our approach as follows [1]: 



N K+ {q,T) = 3x 
N K -(q,T) = 3x 



C M f d 3 p 



(M K F K )W J (2vr)3 e —u(T) + \p- q\/T + 1 .[f^llT ' 

e v s/ +1 

C M f d 3 p 1 1 



{M K F K )W J (2tt)3 e v{T) + \p- q\/T + 1 J^^jt ^ 



N /^ T x o C M f d 3 p 1 1 r , „x 

n±[q ' ] ^ (M^)^J(2nr e -v(T) + \p-q\/T + le v(T)+p/T + 1 > 

where p is a relative momentum of quarks and anti-quarks coalesced in to a meson with 
a 3-momentum q at a temperature T. Since the main contribution to the integrals comes 
from the relative momenta of order p ~ T, so that quark and anti-quarks coalesce at 
relative momenta of order p ~ T. This agrees with the order of transversal momenta 
of hadrons, q± ~ 2 4- 3 T, coupled in the center of mass frame of heavy-ion collisions. 
The factor 3 corresponds the number of quark color degrees of freedom, M# = 500 MeV, 
F K = 160 MeV, M n = 140 MeV and F^ = 131 MeV are the masses and the leptonic 
coupling constants of the K and 7r mesons, respectively [10]. The dimensionless parameter 
Cm is a free parameter of the approach. It is the same for all low-lying mesons. 

The multiplicities of the vector meson production, for example, such as and p 1 ^ 
we define as 



N K * + {q,T) = 3x 
N K *-(q,T) = 3x 



C M f d 3 p 1 1 



{Mk.F k )W J (2vr)3 e -„(T) + \p- q\/T + 1 Jp^jT ^ 
C M [ d 3 p 1 1 



(Mk.F k )W J (2tt)3 e v(T) + \p - q\/T + x Jp^jT ^ ' 

N /^ T x o C M f d 3 p 1 1 n q\ 

VW ' ] (M p F^J (2vr)3 e -u(T) + \p-q\/T + le u(T)+p/T + 1 A ^ ) 

where Mk* = 892 MeV and M p = 770 MeV are the masses of the K* and p mesons, 
respectively [10]. 

lr The numbers of quarks (q, s) and anti-quarks (q,s) and the coefficients of proportionality are free 
parameters of a simple coalescence model. Therefore, a simple coalescence model contains seven free 
parameters. Five of them can be fixed from experimental data [2] . 
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In the case of baryons and antibaryons we suggest to define the multiplicities by using 
the diquark-quark picture of baryons and antibaryons. For example, the multiplicities of 
the proton (p) and antiproton (p) we write in the form 



N p (q,T) = 



3! C 

x 



x 



3! (M p F n ) 3 / 2 
d 3 p 1 



(2tt)3 ( e -v(T) +p/T + x y e -viT) + \p - q\/T + x ' 



3! (M P F W )3/ 2 
y d 3 p 1 1 

X 7 (2^F ( e !/(T) + plT ~ ^ e i/(T) + \p- q\/T + J ' (L9) 

where a momentum p has a meaning of a relative momentum of three-quark (three-anti- 
quark) system, M p = 940 MeV is the mass of the proton and antiproton. As well as in the 
meson case the main contribution to the momentum integrals comes from the momenta 
of order p ~ T providing a coalescence of three quarks (three anti-quarks) into baryons 
(anti-baryons) at the momenta of order p ~ T. That is again of order of transversal 
momenta of the produced hadrons, q±_ ~ (2 3) T, coupled in the center of mass frame of 
heavy-ion collisions. The factor 3! in the numerator is related to the quark colour degrees 
of freedom and defined by SijkS 1 ^ = 3!, where i,j and k are colour indices and run over 
% — 1,2,3 each. In turn, in the denominator the factor 3! takes into account the identity 
of three light quarks (qqq) and three antiquarks (qqq)- In the isotopical limit we do not 
distinguish u and d quarks as well as u and d antiquarks. The dimensionless parameters 
Cb and C B are free parameters of the approach. Each of them is equal for all components 
of octets of baryons and antibaryons, respectively, but Cb ^ C§- 

The paper is organized as follows. In Sect. 2 we calculate the theoretical values of 
multiplicities of the hadron production from the thermalized QGP phase. The theoretical 
predictions and experimental data are adduced in Table 1. In the Conclusion we discuss 
the obtained results. A possible estimate of the absolute values of our input parameters 
is discussed through the application of our approach to the calculation of the number 
of baryons and antibaryons relative to the number of photons at the early stage of the 
evolution of the Universe assuming that this evolution goes through the intermediate 
thermalized QGP phase. 



2 Multiplicities of hadron production from the ther- 
malized QGP phase 

Now let us proceed to the evaluation of multiplicities of hadron production from the 
thermalized QGP phase of QCD. The theoretical predictions for the different ratios of 
hadron multiplicities we compare with experimental data adduced in Table I of Ref. [11]. 
These are the data of various experiments given by NA44, NA49, NA50 and WA97 Collab- 
orations on Pb+Pb collisions at 158 GeV/nucleon. Also we compare our results with the 
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experimental data obtained by NA35 Collaboration on S+S collisions and NA38 Collab- 
oration on O+U and S+U collisions at 200 GeV/nucleon. From Table 1 of this paper one 
can see that in the whole the experimental data for the hadron production are obtained 
for rapidities ranging over the region 2.3 < y < 4.1. The relation between a 3-momentum 
q and a rapidity y reads 



q= \l M 2 sh y + q\ ok y > M shy, (2.1) 

where M and q± are the mass and the transversal momentum of the produced hadron. 
For rapidities y G [2.3,4.1] we get 

q > (5.0 4- 30.2) M. (2.2) 

Thus, for K mesons and hadrons heavier than K mesons typical momenta are of order 
of 2.5 GeV and greater. This gives a possibility to investigate the momentum integrals 
defining multiplicities of the hadron production at q — > oo. As has been shown in Ref. [1] 
the ratios of the multiplicities Rk+k-(<1, T) = N K +(q,T)/N K -(q,T) and i?^+ 7r +(g,T) = 
N K +(q,T)/N w +(q,T) are smooth functions of q, wobbling slightly around the asymptotic 
values obtained at q — > oo, and describe good the experimental data at T = 175 MeV. 
Below the theoretical results on the multiplicities of the hadron production we would 
compare with experimental data at T = 175 MeV. In this case it is obvious that the typical 
momenta of the momentum integrals defining the multiplicities of hadron production are 
of order p ~ T. Therefore, in our approach the momenta at which quarks and anti-quarks 
coalesces into hadrons should be of order p ~ T. This agrees with the order of transversal 
momenta of hadrons, q± ~ 2 4- 3 T, coupled in the center of mass frame of heavy-ion 
collisions. 

Hence, for rapidities ranging over the region 2.3 < y < 4.1 the typical total 3-momenta 
of hadrons are greater than 1 GeV. Thereby, multiplicities of the hadron production can 
be calculated in the asymptotic regime at q ^> T. 

We would like to accentuate that the relative momenta at which quarks and anti- 
quarks coalesce into hadrons are described effectively by the momentum of integration 
p. 

At q ^> T the multiplicities of the hadron production defined by Eqs. (|1.7|) - (|1.9|) can 
be represented in the following form 

N w+ (q,T) = Nw -(q,T)=NAq,T) = -^^e- ( l/ T UT), 

N K+ {q,T) = N K o(q, T) = ,J C " e~^ T e+^ T ) I K (T), 

Nk- (<f> T) = N R0 (q,T) = (M ^ )3/2 e-*l T e-»V) I R {T), 

N K o(q,T) = \ N K0 (q, T)+ 1 - N ro (q, T) = \ (M ^ )3/2 ^ e+ ^ T) (T), 

N (a T) = sin 2 6 3 ° M e~ ^ T I (T) + cos 2 3 ° M e~^ T I(T) 

NAqT) = 3 ° M e-Q/ T UT) 
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N P (q,T) 
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where the structure functions h{T) (i = n, K, K, . . .) are defined by 
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The numerical values of the integrals are obtained at m s = 135 MeV and T = 175 MeV. 

The theoretical ratios of multiplicities of the hadron production which we compare 
with measured experimentally we define as follows 
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R (aT\ N ^ T) ( M p FA 3/2 I^T) 

^°- 0)(g ' T) " nazt) + nazt)-[m^s) ' 
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The constant F$ = 3.5 F w is related to the leptonic constant of the pseudoscalar meson 
containing only s-quarks, ss [12]. We have estimated F$ in agreement with the exper- 
imental data on the r](550)/n and 0(lO2O)/7r production. For the description of the 
multiplicity of the r?(550) meson production we have taken into account that the low- 
energy meson phenomenology [10,12,13] gives the following quark structure of the 77(550) 
meson: 

^(550) = (qq) sin# + (ss) cosO, (2.6) 

where = $o — $p with $ — 35.264°, the ideal mixing angle, and $p, the octet-singlet 
mixing angle. Recent analysis of the value of the octet-singlet mixing angle carried out 
by Bramon, Escribano and Scadron [13] gives dp = — 16.9 ± 1.7°. For the 0(1020) meson 
we have supposed the ss quark structure [10,12]. 



3 Conclusion 

The theoretical and experimental values of the ratios of hadron production are adduced 
in Table 1. From Table 1 one can see a good agreement between presently available set 
of hadron ratios measured for various experiments given by NA44, NA49, NA50 and 
WA97 Collaborations on Pb+Pb collisions at 158 GeV/nucleon, NA35 Collaboration on 
S+S collisions and NA38 Collaboration on O+U and S+U collisions at 200 GeV/nucleon 
and theoretical predictions for the ratios of multiplicities of hadron production from the 
thermalized QGP phase at a temperature T = 175 MeV. Save the ratio A/p, (A/p) t h = 
2.081 and (A/p) exp — 3 ± 1, the theoretical results agree with the experimental ones with 
accuracy better than 18%. 

In the our approach multiplicities of hadron production are defined by momentum 
integrals on quark (antiquark) distribution functions in accordance with the phenomeno- 
logical quark structure of the hadron. For the analysis of the multiplicities of the baryon 
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and antibaryon production in terms of the quark and antiquark distribution functions we 
have followed the diquark-quark picture for baryons and antibaryons. This has allowed 
to describe multiplicities of the baryon, antibaryon and meson production on the same 
footing. 

For the explanation of experimental data on the hadron production in ultrarelativistic 
heavy-ion collisions we have used three input parameters C b /Cb, Cm/Cb and F$. These 
parameters are related to the spatial volumes of hadronization of the quarks and anti- 
quarks from the thermalized QGP phase. The first two parameters have been fixed from 
the experimental data on the ratios (0/fi) exp = 0.46±0.15, (A/K^) exp = 0.65±0.11. This 
gives Cb/Cb = 0.46 ±0.15 and Cm/Cb = 0.23 ±0.04. In turn, the value of the parameter 
Fs = 3.5 F n = 458.5 MeV is a result of a smooth fit of the ratios of hadrons containing the 
ss and ss components in the quark structure. In the bulk our approach to the hadroniza- 
tion from the QGP phase of QCD has succeeded in describing 21 experimental data on 
ultrarelativistic heavy-ion collisions. 

Unlike other approaches [2-4] the ratio of the Q and Q baryon production is an input 
parameter in our model Cb/Cb- The ratio Cl/Q does not depend in our approach on 
both the momenta of baryons and the temperature. The former is due to the zero-value 
of the s-quark chemical potential, fj, s — // g = 0. As a result the ratio Cl/Q can be only 
fitted in our approach. By fitting the ratio fl/fl from experimental data and applying this 
value to the description of other ratios of the baryon and antibaryon production for the 
thermalized QGP phase we have found a good agreement with experimental data. This 
confirms a self-consistency of the approach. 

In turn, the distinction between the parameters Cb and C B can be related to a well- 
known factor of the baryon-antibaryon asymmetry in the Universe which one could put 
phenomenologically at the early stage of the evolution of the Universe whether the baryon 
synthesis in it goes through the intermediate QGP phase. Indeed, as has been stated 
by Borner [29]: Within the standard big-bang model, however, there seems to be little 
chance of achieving a physical separation of baryon and antibaryon phases in an initially 
baryon-symmetric cosmological model. If the baryon number was exactly conserved - as 
it assumed to be the standard model - the small asymmetry necessary for our existence 
must be postulated initially. Grand unified theories offer the possibility of creating this 
small asymmetry from physical processes. 

In our approach the baryon-antibaryon asymmetry at the hadronic level can be re- 
alized phenomenologically in terms of a different rate of hadronization of baryons and 
antibaryons caused by the input parameter C b /Cb = 0.46 ± 0.15 fixed through the ex- 
perimental data on the ratio Cl/Q production in ultrarelativistic heavy-ion collisions. For 
the total number of antibaryons relative to the total number of baryons Nb produced 
for the baryon-antibaryon synthesis at the early stage of the evolution of the Universe at 
a temperature T = 175 MeV [30] and gone through the intermediate thermalized QGP 
phase we predict 

N- C~ 

— - = 0.41 x = 0.19 ±0.06. (3.1) 
Nb Cb 

This result can be supported by a trivial estimate in the approximation of the equilibrium 
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baryon and anti-baryon gases. In such an approximation the ratio N^/Nb is defined by 



Nb = -2p B {T)/T = e -6MT)/T = Q ^ {32) 
Nb 

where a chemical potential fi(T) is given by Eq.( |1.4|) and T = 175 MeV. 

Thus, at the early stage of the evolution of the Universe the number of antibaryons 
should be of order of magnitude less compared with the number of baryons, N B ~ 0.2 Nb- 
According to Borner [31] it is more than enough for the existence of the life in the Universe. 
Recall that the standard approach [31] predicts for every 10 9 antibaryons only (10 9 + 1) 
baryons. As has been stated by Borner: It is to that one part in 10 9 excess of ordinary 
matter that we owe our existence] [31]. 

For the early Universe the total number of baryons and antibaryons was roughly equal 
to the number of photons N p h [31]: 

N B + N B 

^r^ 1 - (3 - 3) 

Since the density of photons is equal to [32] 

^ = ^, (3.4) 

V 7T Z 

where V is the volume of the early Universe, and the density of the total number of 
baryons and antibaryons iVg + Nb calculated in our approach at T = 175 MeV amounts 
to 

Nb + Nb = C B x 0.442 x T 3 , (3.5) 
we can estimate the numerical value of the parameter C B : 

C b ^ 0.55 ±0.08. (3.6) 
This gives the estimate of other input parameters Cg and Cm'- 

C b ^ 0.25 ±0.08, 

C M ^ 0.13 ±0.03. (3.7) 

The analysis of the influence of the input parameter Cb/Cb = 0.46±0.15 on the evolution 
of the baryon-antibaryon asymmetry in the Universe from the early Universe up to the 
present epoch and the formation of a dark and strange matter in the Universe [10] we are 
planning to carry out in our forthcoming publications. 

Now we would like to discuss in more details our approach when compared with 
a simple coalescence model [2]. The main distinction of our approach from a simple 
coalescence model is in the correlation between quarks and anti-quarks coalescing into 
hadrons. In fact, in a simple coalescence model quarks and anti-quarks are uncorrelated 
[2]. This has allowed to introduce separately the number of light quarks q and light anti- 
quarks q and the number of strange quarks s and strange anti-quarks s [2]. Moreover, 
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this has turned out to be of use in order to hide the dependence of quark and anti- 
quark distribution functions on a temperature T and a chemical potential fi(T) in the 
number of light quarks q and light anti-quarks q. A non-vanishing chemical potential of 
strange quarks H S {T) is hidden in the number of strange quarks s and strange anti-quarks 
s. Then, the calculation of multiplicities of hadrons produced from the QGP phase in a 
simple coalescence model resembles a quark counting. In fact, multiplicities of hadrons are 
proportional to products of the number of quarks (q, s) and anti-quarks (q, s) in accord the 
naive quark structure of hadrons. Since quarks and anti-quarks do not correlate, so that 
the multiplicities of hadrons turn out to be independent on the momenta of hadrons. Then, 
the numbers of quarks (q, s) and anti-quarks (q, s) and the coefficients of proportionality, 
the coalescence coefficients C p , C\, C=, Cq and Cp, C^, Cg, Cq, are free parameters 
of a simple coalescence model. Therefore, a total number free parameters appearing in 
a simple coalescence model for the description of baryon and anti-baryon production 
is equal to twelve. By the assumption C p /Cp = C\/C A = Cz/C~ = Cq/Cq = 1 the 
number free parameters has been reduced up to five (q, s, q, s, C), where C is a common 
for baryons and anti-baryons coalescence coefficient. Two of these free parameters have 
been fixed from the fit of experimental data on the baryon and anti-baryon production: 
qjq = 0.41 ±0.02 and s/s = 0.75 ±0.06 [2]. Thus, there are three free parameters left in a 
simple coalescence model applied to the description of baryon and anti-baryon production 
from the QGP phase. It is also important to note that a simple coalescence model [2] 
explains only multiplicities of baryon and anti-baryon production. In fact, save the ratio 
of multiplicities of the K + and K~ mesons none other multiplicities of pseudoscalar and 
vector mesons have been predicted within a simple coalescence model [2]. Therefore, it is 
not completely clear how many free parameters would be added in a simple coalescence 
model for description of multiplicities of pseudoscalar and vector meson production. 

In our approach, where quarks and anti-quarks coalescing into hadrons are correlated, 
we have six free parameters T, /i(T), Cm, Cb, C& an d F s . Five of these parameters 
T = 175 MeV, fi(T) given by Eq.Q)> F s = 3.5 F w = 458.5 MeV, C M /C B = 0.23 ± 0.04 
and Cg/Cs = 0.46 ± 0.15 are fixed. Therefore, only there is one free parameter left in 
the approach. Thus, if to take into account that within our approach we have described 
not only multiplicities of baryon and anti-baryon production from the QGP but also 
multiplicities of pseudoscalar and vector meson production, all together 21 experimental 
ratios, our approach to the thermalized QGP looks much more preferable with respect to a 
simple coalescence model. If to remind that in our approach due to correlations between 
quarks and anti-quarks we are able to follow a dependence of multiplicities of hadron 
production on the hadronic momenta, so an advantage of our approach with respect to a 
simple coalescence model becomes obvious. 
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Table 1. The theoretical ratios of multiplicities of the hadron production are compared 
with the experimental data obtained by NA44, NA49, NA50 and WA97 Collaborations 
on Pb+Pb collisions at 158 GeV/nucleon, NA35 Collaboration on S+S collisions and 
NA38 Collaboration on O+U and S+U collisions at 200 GeV/nucleon. The theoretical 
multiplicities are calculated at the temperature T = 175 MeV. 
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